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Abstract 

We study a variant of bootstrap percolation in which growth is 
restricted to a single active cluster. Initially there is a single ac- 
tive site at the origin, while other sites of Z 2 are independently oc- 
cupied with small probability p, otherwise empty. Subsequently, an 
empty site becomes active by contact with two or more active neigh- 
bors, and an occupied site becomes active if it has an active site 
within distance 2. We prove that the entire lattice becomes active 
with probability exp[a(p)/p], where a(p) is between —ir 2 /9 + Cy/p and 
— 7r 2 /9 + C v ^j(logp -1 ) 3 . This corrects previous numerical predictions 
for the scaling of the correction term. 

1 Introduction 

Local bootstrap percolation is a 3-state cellular automaton on the square 
lattice Z 2 denned as follows. At each time step t = 0, 1, 2 . . ., each site of 
Z 2 is either empty (o), occupied (•), or active (*•). Let p G (0,1). The 
initial configuration is given by a random element a of {o,«,*} z2 under a 
probability measure P p in which 

F p [a(0)=*]=p, P p [a(0) = o] = 1 - p; 

P p [a{x) = m]=p, F p [*(x) = o]=l-p, (x ^ 0); 
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and the states of different sites are independent. (Here = (0, 0) G Z 2 is the 
origin). Subsequently, the configuration at time t + 1 is obtained from that 
at t according to the following deterministic rules. 

(LI) Each • becomes * if it has at least one * within ^-distance 2. 
(L2) Each o becomes * if it has at least two *'s within ^-distance 1. 
(L3) All other states remain unchanged. 

We are interested in indefinite growth, i.e., the event that every site in Z 2 
eventually becomes active. The following is our main result. 

Theorem 1 (growth probability) . There exists constants c, C, po > such 
that for p < po, 

-2A + cyp . , w ^X + Cy/pQogp- 1 ) 3 

exp - — < F„ (indefinite growth) < exp , 

p p 

where A = 7r 2 /18. 



Bootstrap percolation. As suggested by its name, the process we have 
introduced is a close relative of (standard) bootstrap percolation. This 
well-studied model may be described in our setting as follows. Initially, each 
site in the square {1, . . . , L} 2 is independently active with probability p and 
empty otherwise, and all sites outside the square are empty. The configu- 
ration then evolves according to rules (L2,L3) ((LI) being irrelevant), and 
one is interested in the probability I(L,p) that the entire square eventually 
becomes active. 

Bootstrap percolation has been studied both rigorously and numerically 
P, El [3], [I7j, both in its own right and as a tool to analyze other models 
[3 El [7J [H [9] . An important property is that, asymptotically as p — > 0, it 
undergoes a sharp metastability transition as the parameter plogL crosses 
the threshold A. More precisely, as proved in [12], I(L,p) converges to 
or 1 if respectively p log L < A — e or p log L > A + e (also see [T21 [TJ] for 
similar results on related models). However, simulations for moderate values 
of p show surprising discrepancies with this rigorous result [HI El [101 [TB] . 
and our motivation in this article is to advance understanding of this latter 
phenomenon. 

Local bootstrap percolation has been an implicit ingredient in many ar- 
guments involving bootstrap percolation, including those in [21 [101 12], and 
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a variant form appears explicitly in 0, [7j (see the discussion below). The 
proof of the main result from [T2] already implies that for the local model, 
the quantity 

— p log P p (indefinite growth) 

converges to 2A as p — > 0, and the main contribution of this paper is to 
identify the speed of this convergence up to logarithmic factors. 

Heuristically, the dominant mechanism for active sites to take over space 
in bootstrap percolation is the presence of widely separated, and essentially 
independent, local bootstrap percolations. Indeed, for small p, bootstrap 
percolation immediately fixates (stops changing) on the overwhelming pro- 
portion of the lattice; the rare nuclei that facilitate growth hence encounter 
configurations of fixated active sites, which should be not very different from 
occupied sites in the local version. This suggests that I(L,p) makes the tran- 
sition from near to near 1 as L 2 P p (indefinite growth) changes from small to 
large, and thus when 2plogL m —p log P p (indefinite growth). In this subject, 
claims not backed by rigorous arguments have a questionable track record; 
nevertheless, motivated by Theorem [TJ we venture the following conjecture. 

Conjecture 1. In the standard bootstrap percolation model, let L — > oo and 
p — > simultaneously. Then, with A = 7r 2 /18 as above and any e > 0, we 
have 

A - p l / 2 ~ e 

if L < exp then IyL.p) — > 0. 

p 

The complementary bound, namely, 

A - cp 1 ' 2 

if L > exp then I(L,p) — > 1, 

p 

for a small enough c > 0, was proved in [10], so Conjecture [1] states that the 
power of p in the correction term is exactly 1/2. 

Rectangle process. It is natural to consider the following variant of the 
growth model defined by (L1-L3), in which we update states in a different 
order. A rectangle is a set of sites of the form R = {a, c} x {&,..., d} C 
Z 2 . The rectangle process is a (random) sequence of rectangles po Q Pi Q 
p 2 C • • • defined in terms of the initial configuration a as follows. If er(0) = o 
then we set pi = for all i. If cr(0) = * then we set po = {0}, and then 
proceed inductively as follows. If = {a, . . . , c} x {b, . . . , d}, then consider 
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the configuration in which every site in pi is active, every site outside the 
enlarged rectangle pf := {a — 2, . . . , c + 2} x {b — 2, . . . , d + 2} is empty, 
and all sites in pf \ pi have the same state as in a. Apply the update rules 
(L1-L3) to this configuration until the configuration stops changing. It is 
readily seen that the resulting set of active sites is again a rectangle; let p i+ i 
be this set. 

It is straightforward to check that the set of eventually active sites in 
the local bootstrap percolation model is identical to the limiting rectangle 
lim^oopj, and in particular indefinite growth occurs if and only if the latter 
equals Z 2 . 

The rectangle process can be described via a countable-state-space Mar- 
kov chain, whose state represents the current rectangle together with infor- 
mation about which sites on its sides have been examined. In principle, this 
allows for a computational approach to estimating the probability of indefi- 
nite growth; as we do not presently pursue this, we will not give a detailed 
description of the chain. 

Variant models. In addition to the standard model, our methods adapt 
to the following two variants. In the modified local model, we replace rule 
(LI) with: 

(L1F) Each • becomes * if it has at least one * within £°°-distance 1. 

In the Frobose local model, we replace (LI) with: 

(L1M) Each • becomes * if it has at least one * within ^-distance 1. 

These models may be regarded as local versions of modified bootstrap 
percolation [T2l [T3| [15] (in which an empty site becomes active if it has at 
least one active horizontal neighbor and at least one active vertical neigh- 
bour), and Probose bootstrap percolation [9] (in which activation of a 
site requires a horizontal neighbor, a vertical neighbor, and the diagonal 
neighbor between them all to be active). 

Theorem 2 (growth probability for variant models) . Theorem U\ holds for 
the modified and Frobose local models, but with A = ir 2 /G. For the Frobose 
model the upper bound can be improved to exp[(— 2A + C \fp (logp -1 ) 2 ) / 'p\. 

The proof of Theorem [2] follows the same steps as that of Theorem [1] (with 
a few simplifications). We therefore omit the details and instead summarize 
the differences in Section HJ 
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Comparison with numerical results. For the modified and Frobose 
models, the resulting rectangle process, together with the associated Markov 
chain, is much simpler, which makes a computational approach more inviting. 
Indeed, in [7] , using computer calculations together with heuristic arguments, 
the authors obtained the estimates 

-p logP p (indefinite growth) « 2A - 6.22 p a333 

for the modified local model, and 

—p logP p (indefinite growth) « 2A - 5.25 p°' 388 

for the Frobose local model, as p — > 0. These estimates seems to fit well 
down to p ~ 0.01, but Theorem [2] contradicts them in the limit, because the 
estimated powers 0.333 and 0.388 are less than 1/2. 

In a similar vein, in [16] , using interpolation between simulations and the 
rigorous result of [12], it was estimated for the standard (non-local) bootstrap 
percolation model, that the metastability transition occurs at 

phgL^ A- 0.51 p ' 2 . 

This would again be inconsistent in the limit p — ► with our Conjecture [TJ 
These deceptively simple models thus present yet another puzzle of the 
type "why are we not able to see the asymptotic behavior in simulations?" 
That computer simulations can be so misleading [121 [UJ is arguably the 
primary lesson learned from more than two decades of research into boot- 
strap percolation by mathematicians and physicists, and the present article 
is another contribution to this theme. 

Outline of proof. We conclude the introduction with a brief outline of 
the rest of the paper. The bulk of the work, contained in Section [2j is in 
proving the upper bound in Theorem [TJ We consider the rectangle process 
as the side length grows from 1/^/p to p~ l logp -1 . We sample this process 
in a sequence of coarse-grained steps, in which the side length increases by 
a factor of roughly 1 + yfp at each step. These step sizes are chosen to 
balance the entropy factors (related to the number of possibilities) , and errors 
arising from the replacing true probabilities of growth with larger quantities. 
The latter are probabilities of "no double gap" conditions for growth from 
a smaller to a larger rectangle, as in [T2]. One complication, contributing 
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to the logarithmic error factor between our upper and lower bounds, is that 
corner sites between the two rectangles may allow growth in two directions. 

The rectangle process can grow via rectangles of unequal dimensions; 
however, by the variational principle of [12], the most likely scenario is the 
most symmetric one, i.e., growth of squares. This yields the value of the 
constant A. 

The lower bound in Theorem [JJ is a straightforward consequence of meth- 
ods from [10], which we summarize in Section [3j The idea is to consider 
paths in the space of rectangle sizes that deviate from symmetric growth 
by approximately at scale 1/p (and are therefore not prohibitively im- 

probable). The increased entropy from such paths is enough to introduce 
an additional multiplicative factor of exp(c/ v /p) to the bound on the growth 
probability. 

2 Upper Bound 

It will be convenient to assume throughout that the probability p is suffi- 
ciently small, e.g., p < 0.1 suffices. Then we denote 



Note that q = p+p 2 /2+p 3 /3+- ■ ■ for small p. Throughout, we use c, Ci, c 2 , . . . 
and C, Ci, C 2 , . . . to denote positive (respectively small and large) absolute 
constants, which may in principle be explicitly computed. For a rectangle 
R = {a, • • • , c} x {b, . . . , d} we denote its dimensions dim(i?) = (c — a + 
l,d-b + l). 

We will bound the probability of indefinite growth by summing over pos- 
sible trajectories for the rectangle process. However, the total number of 
trajectories is too large, so we need the following concept. Take a sequence 
of rectangles Ri,...,R n+ i, where n > 1. Write dim(_Rj) = (a*, foj) and 
Si := a i+ i — cii and ti := b i+ i — b^. We call the sequence good if it has 
the following properties. 
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-log(l -p); 
[q^logq- 1 ]. 



6 



s 7 


Sq 


s 5 


Ss 


R 


s. 


Si 


s 2 


S3 



Figure 1: The rectangles Si, . . . , S$. 

(i) E Ri C R 2 C ••• C i? n+1 . 

(ii) min(aiA) e [A, A + 3]. 

(iii) a n + b n < B. 

(iv) a n+ i + > B. 

(v) For f = 1, . . . , n, either Sj > a iy /q or U > biy/q. 

(vi) For i — 1, . . . ,n, both Si < ai^/q + 4 and ^ < b iy /q + 4. 

We next define some useful events. The columns of a rectangle R = 
{a, . . . ,c} x {b, . . . ,d} are the sets {a} x {b, . . . , d}, . . . , {c} x {b, . . . , d}. We 
say that R has a double gap in the columns if two consecutive columns 
are entirely empty in the initial configuration a. Double gaps in the rows 
are defined similarly. For a rectangle R, define the event 

G(R) := {R has no double gaps in the columns or rows}. 

Moreover, we recall the following definition from [12]. For two rectangles 
i? C ]?', define subrect angles Si, . . . , S 8 (some of which may be empty) ac- 
cording to Figure [fl so that R' is the disjoint union of R, Si, . . . , S 8 . Define 
D(R,R') to be the event that each of the two rectangles Si U S$ U S 7 and 
63 U S4 U S5 has no double gaps in the columns, and each of the two rectangles 
Si U S 2 U S 3 and £7 U S e U S5 has no double gaps in the rows. The idea is 
that this event is necessary for the growth to proceed from R to R' . 

Finally, define the event 

R contains and has one dimension in 1 
[B - A - 10, B - A] and the other in [1, A] J ' 
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Lemma 3 (key inclusion) . If indefinite growth occurs, then either E occurs, 
or else there exists a good sequence Ri, . . . , R n +i of rectangles such that 

n 

G(R 1 )nf]D(R i ,R i+1 ) 

occurs. 

Proof. We start with the following observation: if the rectangle process ever 
encounters a given rectangle R, then G(R) occurs. To see this, suppose on the 
contrary that R has two consecutive vacant columns, say If these columns 
contain the origin then no growth occurs. Otherwise, so long as the growing 
rectangle includes no site above or below R, it contains no site in the two 
columns. This follows by induction on the steps of the rectangle process (no 

in the two columns can have two adjacent *'s within R, while no • outside 
R can contribute to a site becoming * without first becoming * itself). 

Now define the good region for dimensions of rectangles to be 

T = {(a,b) £ Z 2 : a,b > A and a + b < B) 

(see Figure [2]). Note that in the sequence of dimensions (dim(pj))j>o for the 
rectangle process, each coordinate increases by at most 4 at each step, while 
if indefinite growth occurs then at least one coordinate increases by at least 

1 at every step, and both coordinates tend to oo. 

Therefore, assuming indefinite growth, if (dim(pj))j>o never enters T, 
then E must occur, since the sequence must "escape" near a corner of the 
good region. On the other hand, if the sequence does enter T, then let R\ be 
the first pj such that dim(p,,) £ T. Then property (ii) of a good sequence will 
be satisfied. Thereafter, define R2, . . . , R n +i iteratively as follows. Assuming 
dim(.Ri) = (dj, bi), let R4+1 be the first rectangle encountered by the rectangle 
process after Ri such that either Sj > a iy /q or U > b iy /q (where (si,U) = 
dim(i?j + i) — dim(i?j)). This ensures (v),(vi) are satisfied. Continue in this 
way, stopping at the first such rectangle, R n +i, whose dimensions are outside 
T. Then (iii),(iv) are satisfied. 

The rectangles thus constructed form a good sequence. From the above 
observations, fl" = i 1 1 G(-Ri) occurs, and since D(R,R') C G(R') the result fol- 
lows. □ 
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[B - A, A) 



Figure 2: The sequence of dimensions (a, b) of the rectangle process (bold 
curve), as it passes through the good region T. The small rectangular window 
has dimensions (a iy /q,biy/q). 



Next we bound the probabilities of the events appearing in Lemma [31 We 
introduce the following functions from [T2] : 



u+ y/u(4:-3u) 
P( u ) = -~2 and 9{z) = -log/3(l -e ). 

The threshold A is determined by the integral (see [T2], and also [13] ) 



/•oo ^2 

J 3{z)dz = X = -. 



(1) 



Lemma 4 (double gaps). If R is a rectangle with dim(i?) = (a,b), 

f p (R has no double gaps in the columns) < exp [ — (a — l)g(bq)]. 
For rows we have the same bound with a and b exchanged. 
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Proof. See [I2J Lemma 8(i)]. 



□ 



The following is an enhanced version of [T2l Proposition 21]. The reader 
may want to consult that proof for more details. 

Lemma 5 (border event). For rectangles R C R' of dimensions (a,b) and 
(a + s,b + t), 

F p [D(R,R')] 

< exp ( - sg{bq) - tg{aq) + 2[g{bq) + g{aq)\ + stq e ^(b9)+g(m)]\ < 

We remark that the second appearance of the constant 2 is crucial; the 
lemma would not hold with a smaller constant, while a larger constant would 
not suffice for the remainder of our argument, owing to the behaviour of g 
near zero (Lemma [TU1 below) . 

Proof of Lemma We first claim that 

F P [D(R,&)] < (fc V(* - p) st - k e-^- 2k - 2 ^e^ t - 2k - 2 ^ aq \ 

k=0 ^ ' 

To prove this, we split the event according to the set of occupied sites in the 
set Si U S 3 U S5 U S 7 . If the number of such sites is k, they divide S 4 and 
S 8 into at most k + 2 contiguous vertical strips, each of which has no double 
gaps in the columns. A similar argument applies to rows in S2 and Sq, and 
the relevant events are independent. Then use Lemma HI 

From the above, dropping the power of (1 — p) and using the binomial 
expansion, we obtain 

Fp[D(R,R')] < e - s 9{bq)-tg{aq)+2[g{bq)+g{aq)] ^ _|_ p e 2[g(bq)+g(aq)] ^ st 

Applying 1 + z < e z and p < q to the last factor finishes the proof. □ 

Lemma 6 (escape probability). For some absolute constant c > and all 
p < 0.1, 

F P (E) <exp [-cg^Qogg- 1 ) 2 ] . 
Proof. By Lemma H] and the definition of E, 

F p (E) < \\A 2 B ■ exp[(B -A - ll)g(Aq)}. 
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The factor in front of the exponential comes from the number of possible 
choices of rectangles. Now use the definitions of A and B and the fact that 



g(e) ~ - \oge 



-i 



as e 



0. 



□ 



We have provided all the probabilistic bounds we need, and we now pro- 
ceed to analytic estimates of the expressions that appear. The following will 
be applied to the first two terms in the exponential in Lemma [H 

Lemma 7 (variational principle). Let A,B be any integers with A > 4 and 
B > 2A, and suppose (aj, 6j)j = i v .. jn+ i satisfy = b = A, and Sj := a i+1 — 
ai > and t,i := b i+ \ — bi > 0, and properties (ii)-(iv) in the definition of a 
good sequence. For any q > and any positive, smooth, convex, decreasing 
function g : (0, oo) — > (0, oo), 



Proof. We make use of the following technology from [121 Section 6]. If 7 is 
a piecewise-linear path in the quadrant (0, oo) 2 , with each segment oriented 
in the non-negative direction of both co-ordinates, we define 



Let 71 be a piecewise-linear path with vertices (a q, b q), . . . , (a n+ iq, b n+ iq) 
(in that order). Since g is decreasing, the sum in the statement of the lemma 
is at least w{^i)/q (compare [12], Proposition 16]). 

Also let 72 be the straight line path from (a n+ iq, b n+ iq) to (Bq,Bq), 
and let A be the straight path from (Aq, Aq) to (Bq,Bq). The variational 
principle in [121 Lemma 16] states that w is minimized by paths that follow 
the main diagonal, therefore 




i=0 
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w(ji) + u>(7 2 ) > w(A). 




On the other hand, 

w( l2 )<2Bq-g(Bq/2), 
because min{a n+ i, b n+ i} > B/2 and (B — a n+ i) + (B — b n+ i) < 2B. 



□ 
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Lemma 8 (integral). With q,A,B,g as defined previously, forp < 0.1 and 
an absolute constant C\, 



f 9 7T 

/ g >—- C ly /q log q' 1 . 
JAq lo 

Proof. Use (Tj[|) together with the asymptotics 



/ € 1 , _i , f°° K-^oo 1 _2R- ■-, 

# ~ 2 e ge y ^ ~ 2 e 

Lemma 9. For some C*2 and aZ/p < 0.1 we have Bg(Bq/2) < C^logg -1 . 
Proof. As g(if ) ~ e~ 2 ^ as if — > oo, and by the definition of B, 

Bg{Bq/2) < ^ logg"^' • e" 2 '! log ^ 1 " 1) . □ 

The next two lemmas will be used to bound the last two terms in the 
exponential in Lemma El 

Lemma 10. For every a < B we have e 2g( - a ^ < ^ logo -1 . 

Proof. We have e 9 ^ = l//3(l - e~ z ). Moreover, (3{l - e~ z ) ~ z 1 ' 2 as z -> 
and j3(l — e~ z ) — > 1 as z — > oo. It follows that, for a large enough M, 



SU P tttt^ V < 2vM. 

o<z<m p(l - e~ z ) 

Therefore, 

Lemma 1 1 (summation bound) . Let n and , 6j (z = 1 , . . . , n + 1 ) be positive 
integers and suppose that Sj := a i+1 — aj > and ij := — 6j > for i = 
1, . . . ,n. Further, assume that (ii)-(vi) in the definition of a good sequence 
are satisfied. Then 

1 n t 

n < ^-Xogq' 1 and } < C 4v /glogg _1 . 
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Proof. To bound n, use (v),(ii) and (iii) to get 

Then use this to prove the second bound as follows. As dj > A and bi > A 
and A^fq > 1, (vi) implies Sj < 5a iA /g and tj < 5bi^/q. Therefore 

n « v- 

y ^ < n . 25g. □ 

i=i 

It will be convenient to bound the probability of G(R\) (associated with 
the entry of the rectangle process into the good region) in terms of a factor 
that can be easily combined with those obtained in Lemma [5], in order to use 
Lemma [3 The following lemma deals with this estimate. 

Lemma 12 (entry term). Let R\, . . . , R n +i be a good sequence of rectangles 
and let = bo := A, with s := a± — a and t := bi — 6 - Then, for some 
absolute constant C^, 

P P [G(Ri)] < exp [ - s g(b q) - t g(a q) + C^ 2 log q' 1 ] . 

Proof. Without loss of generality suppose a\>b\. By Lemma IH 

P[G{Rx)\ < e^ 01 " 1 ^*). 

Now bi < A + 3. We have g'(z) ~ l/(2y/z) as z — > 0, and g is convex. 
Therefore 

< g(Aq) - g(b iq ) < 3q ■ g'(Aq) < C^. 

Also to < 3, so 

e tog{aoq) < e Zg{Aq) < C 7 g~ 3 / 2 . 

Therefore, 



exp [-s g(b q) - t g(a q)} 
Lemma 13 (entropy). The number of good sequences of rectangles is at most 

'a 



exp 



^(logg- 1 ) 2 
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Proof. There are at most B 4 choices for R\. For each subsequent rectangle 
there are at most (5y / g + 4) 4 choices, corresponding to the distance by which 
the rectangle grows in each of the four directions. The number of steps is 
n+ 1, which is bounded by Lemma [TT1 This gives the following upper bound: 



4(q- 1 / 2 logg- 1 +l) 



S 4 (IV^ + 4) 



We are now ready to prove the upper bound in the main theorem. 



□ 



Proof of Theorem^ (upper bound). Take p < 0.1. For a good sequence of 
rectangles, the events G(R\), D(Ri, Ri+i) appearing in Lemma [3] are all in- 
dependent, since they involve the initial states of disjoint sets of sites. There- 
fore, 



P p (indefinite growth) < F P (E) + F p i G ( R i)} H P p [ D ( R ^ ^+0] > 



Rl, — ,Rn+l 



i=l 



(2) 

where the sum is over all possible good sequences of rectangles. We now 
proceed to bound the various terms. 

We start by using Lemmas 151 and [TOl to bound F p [D(Ri, Ri+i)] above by 



exp[-Sig(biq) -Ugfaq)] f^logg M exp 

C 2 

< exp [-sigipiq) - Ug(aiq)] ■ (log g" 1 ) • exp 



SiUq — logg — logg 



aiq kq 

Cg / -l\2 S{ti 

(logg J 



dibi 



Using this and Lemma [U we next bound F p [G(Ri)]Yl" =1 Wp[D(Ri,Ri + i)] 
above by 



exp ( ~ S [ 8 *9{ b iQ) + Ug{aiq)] J 



x exp 



9l 



logg 1 



— (logg ) 



■ exp 



^9 /, _1\2 v— v SiU 
— (log q ) > — 



i=l 



a-ik 



By Lemma dU this expression is at most 

exp ( - y^\sig{biq) + Ugfeq)] ) • exp 

14 



i=0 



1\3 



ClO /, 



Next, we use Lemmas [7J M and [HI to bound 



n 2\ (J 
Y] \sig{hq) +tig(aiq)] > JMogg -1 . 

Finally, we substitute the last two bounds, together with Lemmas [13] and 
El into ([2D to get 



Pp(indefinite growth) < exp [—cq 1 (logg x ) 2 ] 



(logg 



-1\2 



■ exp 



2A d 



Mogg 1 



exp 



C 
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(log? 



-1\3 



+ exp 
Therefore 

P p (indefinite growth) < exp 

L y 

As a final step, use the asymptotics q = p + 0{p 2 ) as p —> to replace g with 



— + — logg^ 1 3 
9 



□ 



3 Lower Bound 

Proof of Theorem^ (lower bound) . This follows easily from the proof of 
Proposition 8 in [10]. In that article, it is proved that there is a certain 
event J 7 , satisfying 

P P (^)>exp[-2A/p + c/Vp], 

for small p, and defined in terms of the initial states in the quadrant Q : = 
{0, 1, . . .} 2 C Z 2 , on which the entire quadrant Q eventually becomes active 
in the bootstrap percolation model. On the same event (but replacing * with 
o for sites i^O), the entire quadrant becomes active in the local bootstrap 
percolation model also. This fact would not follow for an arbitrary event, but 
it holds for the particular event T because it is tailored to produce growth 
starting from the origin. 

For the reader's convenience we briefly summarize the construction of the 
event J 7 , referring to [10] for the details. We consider two ways to grow from 
active square [0,a] 2 to the active square [0,6] 2 . The first, symmetrical, way 
is that the top and right sides of the growing square always each have an 
occupied site within distance 2. The alternative, deviant, possibility is that 
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vertical growth is temporarily halted by two adjacent vacant rows, and the 
square first elongates horizontally to the rectangle [0,b] x [0,a], when it can 
finally proceed to grow vertically until it reaches [0, b] 2 . A key estimate, [T0| 
Lemma 13], states roughly that, for a, b on the scale 1/p, the quotient of prob- 
abilities of the second and the first event is at least c\ p exp [— 2C p{b — a) 2 ]. 
This suggests that one should keep b — a of order 1/y/p to incur a deviation 
cost c 2 p. Indeed, one can make about m := c 3 / ^fp deviant steps within these 
constraints, and a careful organization ensures that different choices of these 
steps give disjoint events. The resulting lower bound on P p (^ r ) is greater 
than exp [— 2X/p] by at least the factor 

CO (i/ ^) m • {c2p)m w (c2/csr ' 

which gives the required bound provided we choose c 3 <C c 2 . 

Finally let Q be the event that every semi-infinite line of sites, started 
at arbitrary site and moving in any direction parallel to one of the axes, 
contains at least one • initially. Clearly, P p ((7) = 1 and 

J~ H Q C {indefinite growth}, 

completing the proof. □ 

4 Modified and Frobose models 

In this section we describe the proof of Theorem[2j The proof follows the same 
steps, and in fact a few simplifications are possible. We therefore summarize 
the differences. 

For both the modified and Frobose models, we replace the condition of 
"no double gaps" in the definitions of the events G(R) and D(R, R') with the 
condition that no columns (or rows) are initially entirely empty. The bound 
in Lemma H] then becomes exp[— af(bq)], where 

/(^) = -log(l-e- 2 ). 

The function / then replaces g throughout, and the threshold A then arises 
as J °° / = 7r 2 /6. We also modify the definition of the rectangle process so 
that p + is p enlarged by only 1 in each direction. 
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For the modified model, the change to Lemma H] allows the second 2 in 
the exponential in Lemma [5] to be replaced with 1. In fact this change is 
needed for the argument to go through, because of the different behavior of 
/ near zero. Specifically, Lemma [10] holds with e^ aq ^ on the left side (and 
the same right side). 

For the Frobose model we can make a further simplification to the argu- 
ment, as horizontal and vertical progress now occur disjointly. To be more 
precise, D(R, R') is now defined to be the event that the rectangle process 
makes a transition from R to R' . Then, referring to Figure 1, the two events: 

Si U S2 U S3 and S5 U Sq U S7 have no empty rows, and 
Si U S7 U Sg and S3 U S4 U S5 have no empty columns, 

occur disjointly and are thus negatively correlated by the Van den Berg 
Kesten inequality (see e.g. [H]). This fact eliminates the last summand in 
the exponential in Lemma O makes Lemma HO] unnecessary and reduces the 
power of logp" 1 to 2 in the final bound. 

Open Problems 

(i) Is a power of logp -1 in the upper bound of Theorem [U really necessary? 

(ii) In the modified and Frobose models, can the discrepancy between the 
asymptotic power 1/2 and the numerical estimates, which are closer to 
1/3, be explained? 

(iii) Is it possible to make the same arguments work if the distance in rule 
(LI) is altered to something larger than 2? 
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